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Electronic  attack  on  air  defense  radar  networks  by  a  team  of  Electronic  Combat  Air 
Vehicles  (ECAVs)  is  explored.  The  scope  of  the  attack  considered  is  the  creation  of  a 
coherent  phantom  track  in  two  dimensions — range  and  azimuth — by  a  team  of  coopera¬ 
tively  controlled  ECAVs,  where  each  ECAV  is  capable  of  intercepting  and  sending  delayed 
returns  of  radar  pulses.  Severe  restrictions  are  shown  when  a  constant-speed  ECAV  and 
constant-speed  phantom  track  are  assumed.  Dynamic  limitations  on  an  ECAV  are  pre¬ 
sented  mathematically  and  through  simulation  for  straight  and  circular  phantom  tracks 
and  specified  ranges  on  the  phantom  target  speed,  ECAV  speed,  and  ECAV  region  of  an¬ 
tenna  operation.  Based  on  results  for  a  straight  or  circular  phantom  track,  generalized 
bounds  for  the  initial  conditions  and  time-dependent  flyable  ranges  of  a  team  of  ECAVs 
are  presented.  These  bounds  are  then  used  with  coordination  functions  to  formulate  a 
decentralized  cooperative  control  problem  for  the  scenario  of  three  radars  versus  three 
ECAVs  generating  a  straight  phantom  track. 


Nomenclature 

B  Space  where  placement  of  a  phantom  track  is  possible  for  all  radar/ECAV  pairs 
7g  Set  of  all  phantom  track  functions  f[x)  possible  for  ECAV  team  to  generate 

a  Horizontal  distance  from  center  of  circular  phantom  track  to  radar 

a e  ECAV  acceleration — time  derivative  of  velocity 

b  Vertical  distance  from  center  of  circular  phantom  track  to  radar 

c  Speed  of  light 

circle  2-D  Circular  boundary  of  radius  Rmax  inside  which  radar  can  operate  and  see  targets 
cone  2-D  Wedge-shaped  boundary  inside  which  ECAV  antennas  can  deceive  radar 
E  ECAV 

/  Phantom  track  function  of  x  that  ECAV  team  could  generate  for  radar  network 
h  Coordination  function  providing  individual  ECAV’s  cost  versus  different  phantom  tracks 
He  ECAV  heading  angle — absolute 

hr  Function  that  determines  valid  initial  conditions  and  flyable  ranges  for  ECAV 

hr  Function  that  maps  parameters  to  all  possible  phantom  tracks  for  ECAV  team 

I  Ordered  set  of  index  numbers 

icost  ECAV’s  cost  to  move  to  its  initial  position  required  for  a  given  phantom  track 
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n  Number  of  radars  in  a  given  scenario 
R  Phantom  target  range  from  radar 

r  ECAV  range  from  radar 

Rmax  Maximum  operational  range  of  radar 

rd  Location  of  radar  in  x  and  y  coordinates 

s  Switching  function/indicator 

T  Phantom  target 

t  Time 

td  Time  delay  for  ECAV  to  hold  radar  pulse 
ts  Switching  time 

tcost  ECAV’s  cost  over  time  to  generate  a  given  phantom  track 

x  Horizontal  coordinate  in  2-D  space 

y  Vertical  coordinate  in  2-D  space 

ECAV  Electronic  Combat  Air  Vehicle 

LOS  Line  of  sight  from  radar  to  ECAV/phantom  target 

PRF  Pulse  repetition  frequency 

Subscripts 

0  Value  at  the  initial  time 

c  Variable  set  constant 

i  Index  number  for  radar/EC AV  pairs 

j  Index  number 

to  Nominal  or  average 

psmax  Pseudo-maximum 

psmin  Pseudo-minimum 

x  ^-component  of  variable — horizontal 

y  y-component  of  variable — vertical 

Symbols 

a  Ratio  of  phantom  target  speed  to  ECAV  speed 

(3  Absolute  heading  of  phantom  track — independent  of  radar/EC  AV  pair 

X  Angle  between  current  LOS  and  line  from  center  of  circular  phantom  track  to  radar 

cf)E  ECAV  course — relative  to  9 

<j>T  Phantom  target  course — relative  to  9 

%l>  Angle  between  initial  LOS  and  phantom  target  heading 

p  Radius  of  circular  phantom  track 

9  Constant  angle  from  0  degrees  heading  to  initial  LOS 

9  ECAV/phantom  target  bearing  from  radar 

ve  ECAV  speed 

vt  Phantom  target  speed 

£  Heading  angle  from  center  of  circular  phantom  track  to  radar 
Superscripts 

0  Solved  forward  in  time  using  initial  values 
A  Solution  corresponding  to  ECAV  system  A 
/  Solved  backward  in  time  using  final  values 

I.  Introduction 

This  paper  focuses  on  using  unmanned  Electronic  Combat  Air  Vehicles  (ECAVs)  to  cooperatively  deceive 
a  radar  network  by  creating  a  coherent  phantom  track,  which  allows  the  network  to  detect  and  track  the 
motion  of  an  air  vehicle  that  does  not  actually  exist.  A  radar  network  is  defined  as  a  geometric  arrangement 
(possibly  changing  with  time)  of  two  or  more  radars,  which  are  able  to  communicate  with  each  other;  hence, 
multiple  radars  may  be  used  to  track  and  correlate  the  same  target.  Ref.  1  provides  a  thorough  background  on 
radars  and  how  they  work.  For  this  paper,  generating  a  viable  coherent  phantom  track  requires  a  specified  set 
of  the  networked  radars  to  detect  a  consistent  position  (range  and  azimuth)  and  range  rate  for  the  phantom 
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target  as  it  moves  along.  The  methods  considered  herein  to  accomplish  such  a  task  are  restricted  to  range- 
delay  techniques  applied  through  the  radar  mainlobe.  Obviously,  for  a  team  of  ECAVs — generally  one  per 
radar — attempting  to  deceive  a  radar  network  in  such  a  way,  the  flyable  trajectories  of  each  individual  ECAV 
are  restricted  to  some  exclusive  range. 

Prior  work  by  Pachter  and  Chandler  in  Ref.  2  presents  a  thorough  introduction  to  the  cooperative 
electronic  attack  problem  and  mathematically  describes  the  dynamics  of  a  constant- speed  ECAV  given  a 
straight  or  circular  phantom  track  with  constant  speed.  This  paper  begins  by  showing  through  simulation 
that  assuming  an  ECAV  has  constant  speed  severely  limits  the  range  of  flyable  ECAV  trajectories  given 
a  straight  constant-speed  phantom  track.  In  response,  more  general  theory  for  modeling  ECAV  dynamics 
mathematically  and  developing  bounds  for  ECAV  trajectories  is  contributed,  which  allows  for  ranges  on  the 
ECAV  and  phantom  track  speeds  as  well  as  limitations  on  the  ECAV  antenna  operation  regions.  With  these 
specifications,  the  flyable  ranges  for  an  ECAV  given  a  straight  or  circular  phantom  track  are  developed  and 
verified  via  simulations.  Also,  generalized  mathematical  bounds  for  a  team  of  ECAVs — independent  of  the 
phantom  track  types  considered — are  provided.  One  formulation  of  the  decentralized  cooperative  control 
problem  is  supplied  as  an  example,  which  uses  these  generalized  bounds  on  ECAV  initial  conditions  and 
flyable  ranges  to  determine  the  ECAV  team  cost  versus  different  phantom  tracks  in  a  decentralized  fashion. 

II.  An  Explanation  of  the  Deception  Problem 

ASSUMING  that  an  ECAV  knows  the  maximum  operational  range,  Rmax >  and  location  of  a  radar  with 
pulse-to-pulse  agility  or  discrimination,  the  ECAV  can  intercept  and  appropriately  delay  the  return  of 
the  radar’s  transmitted  pulses  so  that  the  radar  sees  a  phantom  target  beyond  the  ECAV  but  closer  than 
Rmax ■  The  capability  of  digitally  storing,  altering,  and  returning  encoded  pulses  so  that  they  correspond 
to  a  desired  range  and  range  rate  is  formally  known  as  Digital  Radio  Frequency  Memory  (DRFM),  but  will 
often  be  referred  to  simply  as  range  delay  since  it  is  assumed  that  digital  means  will  be  used  to  accomplish 
this  task;  see  Ref.  1  for  more  information.  To  deceive  the  radar  using  only  range  delay  requires  that  the 
ECAV  be  in  the  radar’s  mainlobe;  of  course,  the  phantom  target  must  also  be  on  the  line  of  sight  (LOS) 
from  the  radar  to  the  ECAV.  Thus,  the  phantom  target’s  track  is  a  function  of  the  ECAV’s  range  delay 
and  bearing  from  the  radar.  If  the  radar  operates  at  a  pulse  repetition  frequency  (PRF)  where  pulse-delay 
ranging  is  performed,  then  accurately  delaying  the  radar’s  pulses  is  crucial.  If  the  radar  operates  at  higher 
PRFs  in  track  mode,  then  range  information  is  still  important,  but  the  ECAV  can  focus  on  sending  returns 
with  accurate  Doppler  frequencies  or  range  rates  for  the  phantom  target. 

Keep  in  mind  that  a  phantom  track  flying  roughly 
perpendicular  to  a  radar  will  not  produce  significant 
Doppler  frequencies  and  may  even  get  filtered  out  or 
ignored.  Figure  1  illustrates  how  four  ECAVs  could 
cooperatively  create  a  single  phantom  track  to  deceive 
a  network  of  four  radars  by  using  range-delay  tech¬ 
niques.  For  this  example,  where  all  four  of  the  radars 
share  track  files,  unless  all  radars  see  the  same  phan¬ 
tom  track,  the  track  is  dismissed  as  spurious. 

One  problem  of  interest  is  to  determine  the  allow¬ 
able  trajectories  for  the  ECAVs  given  their  initial  con¬ 
ditions  and  a  time-dependent  phantom  track.  Observ¬ 
ing  figure  1,  one  can  see  that  each  ECAV  behaves  much 
like  a  bead  on  a  string  that  is  rotating  at  some  vari¬ 
able  rate;  the  ECAV  may  slide  up  and  down  freely  but 
must  rotate  with  the  LOS  from  the  radar  to  the  phan¬ 
tom  track.  If  the  speed,  heading,  or  another  variable 
of  the  ECAV  is  limited,  then  it  may  no  longer  slide 
freely;  in  fact,  an  allowable  trajectory  for  the  ECAV 
might  not  exist.  This  idea  can  also  be  expressed  using 
the  concept  of  degrees  of  freedom  (DOF).  Any  general 
phantom  track  in  the  x-y  plane  has  two  DOF,  which  may  be  represented  in  polar  coordinates  by  R  and  6 
with  reference  to  a  radar  as  the  origin.  Likewise,  two  DOF  also  represent  any  free  ECAV  trajectory  in  the 


Figure  1.  Cooperative  mainlobe  deception  of  a  radar 
network  by  creation  of  a  phantom  track. 
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same  plane  and  may  be  represented  by  r  and  9  (see  figure  2  below).  Thus,  an  ECAV  trying  to  generate  a 
given  phantom  track  will  have  one  constraint,  9 ,  and  one  DOF.  This  DOF  is  then  constrained  by  setting  any 
other  ECAV  variable. 

Mainlobe  as  opposed  to  sidelobe  deception  has  the  primary  advantage  in  that  it  is  less  sophisticated  to 
execute.  Its  main  disadvantage  is  the  ability  of  producing  only  one  phantom  track  when  the  number  of 
ECAVs  is  equal  to  the  number  of  networked  radars.  Other  issues  that  need  to  be  addressed  for  mainlobe 
deception  include: 

•  Inaccuracy  of  radar /ECAV  positions  and  time  delay — an  estimation  problem 

•  Limitations  on  ECAV  dynamics  due  to  bounded  ranges  on  its  speed  and  antennas  and  the  phantom 
target  speed 

•  Choosing  the  “best”  phantom  track  for  a  team  of  ECAVs  — a  decentralized  cooperative  control  problem 

•  Processing  power  required  by  DRFM — can  be  better  at  low  PRF 

•  Electronic  requirements  for  generating  returns  with  sufficiently  accurate  range  and  Doppler  frequency 
information. 

The  first  issue  listed  above  is  the  subject  of  current  research  but  is  not  addressed  below.  The  second  and 
third  issues  will  be  addressed  later  in  more  detail. 

III.  Results  of  a  Constant-Speed  ECAV  and  Phantom  Track 

Understanding  first  the  dynamic  limitations  imposed  on  an  ECAV  when  its  speed  and  the  phantom 
track  speed  are  assumed  constant  will  help  motivate  the  later  application  of  realistic  bounded  ranges 
for  the  speed  of  the  ECAV,  its  region  of  antenna  operation,  and  the  speed  of  the  phantom  track.  Using 
mainlobe  deception  as  explained  above,  assume  for  now  that  the  ECAV’s  one  DOF  is  constrained  by  a 
constant  speed.  The  inverse  problem  is  now  of  interest:  given  a  time-dependent  phantom  track  and  an 
ECAV’s  initial  position,  synthesize  the  ECAV  trajectory  required  to  create  the  desired  phantom  track.  For 
ease  of  comparison,  the  following  non-dimensional  variables  are  used,  with  ve  =  1,  Ro  =  1,  and  9q  =  0. 


o  |  ud 

r 

Ro 

VT 

TtO 

ve 

Given  these  definitions  and  a  constant-speed  constant-course  phantom  track  for  this  analysis,  figure  2  below 
illustrates  the  appropriate  variables  and  their  relations.  Definitions  of  the  basic  ECAV  and  phantom  track 
variables  are  also  given  on  the  right  side  of  this  figure. 

The  basic  equations  of  motion  for  the  ECAV  and  phantom  target,  using  polar  coordinates,  are  as  follows. 
Every  derivative  is  taken  with  respect  to  time. 


t  —  cos  <j}E ,  r(0)  =  r0 

(1) 

9  =  -  sin  <j>E,  91 0)  =  0 

r 

(2) 

R  =  acoscj)T,  R{  0)  =  1 

(3) 

9  =  a-sinc/)T,  9(0)  =  0 

(4) 

Without  loss  of  generality,  9  is  assumed  positive.  Manipulating  equations  (1)  and  (2)  results  in  the  following 
differential  equations.  These  same  equations  easily  follow  by  using  the  Pythagorean  Theorem  on  figure  2. 

f  =  yj  1  -  (r9)2,  r(0)  =  r0  (5) 

r  =  -\/l  -  ( r9 )2,  r(0)  =  r0  (6) 
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\)T=  a 


Figure  2.  Variables  and  relations  for  a  constant-speed  ECAV  and  a 
constant-speed  constant-course  phantom  track. 


E  ECAV 

ve  ECAV  speed 

(j>E  ECAV  course 

r  ECAV  range  from  radar 

9  ECAV/phantom  target 

bearing  from  radar 

T  phantom  target 

R  phantom  target  range 
from  radar 

vt  phantom  target  speed 

4>t  phantom  target  course 

ip  phantom  target  heading 


For  both  equations,  the  ECAV’s  range,  including  r0,  must  satisfy  the  following  inequality  with  6  as  a 
function  of  time  for  a  solution  to  exist. 

r(t)  <  i  Vt  (7) 

9 

If  condition  (7)  is  not  satisfied,  the  solutions  to  equations  (5)  and  (6)  become  imaginary.  Physically,  condi¬ 
tion  (7)  makes  it  clear  that  an  ECAV  cannot  fly  at  a  constant  speed  of  one  when  its  range  is  such  that  r6 

is  greater  than  one  (see  figure  2  above).  If  rO  is  less  than  one,  the  ECAV  may  fly  at  a  constant  speed  of 

one  by  having  the  correct  f  component  via  equation  (5)  or  (6).  If  either  equation  (5)  or  (6)  reaches  zero  in 
finite  time  at  say  ts,  the  necessary  condition  for  switching  from  one  equation  to  the  other  at  ts  is  as  follows, 
assuming  the  second  derivative  of  9  exists. 

0(ts)  =  0  (8) 

Finally,  the  equations  defining  the  constant-speed,  constant-course  phantom  track  in  polar  coordinates 
are  as  follows  (see  figure  2  above). 


R(t)  =  \J\  +  a2t2  —  2at  cos  ip 

(9) 

.  .  ( at  sin  ip\ 

9<‘)  =  "CSm(  R(t)  ) 

(10) 

■  a  sin  ip 

1  +  a2t2  —  2 at  cos  ip 

(11) 

To  solve  for  an  ECAV  trajectory  given  this  phantom  track,  equation  (11)  is  inserted  into  equation  (5)  or 
(6),  and  the  solution,  r(t),  is  plotted  versus  9(t ),  which  is  obtained  from  equations  (10)  and  (9).  A  good 
portion  of  the  mathematical  development  given  above  is  from  Ref.  2. 

Figures  3  and  4  below  effectively  explore  equations  (5)  and  (6)  by  plotting  ECAV  trajectory  solutions  for 
various  initial  conditions.  For  each  equation/figure,  the  initial  conditions  are  selected  to  show  the  minimum 
and  maximum  r0  values  yielding  a  ffyable  ECAV  trajectory  for  a  full  90  degrees,  and  also  one  initial  condition 
or  ro  value  where  the  trajectory  runs  into  the  bound  represented  by  condition  (7)  before  9  =  90  degrees 
is  reached.  In  both  figures,  the  maximum  valid  value  for  ro  yields  an  ECAV  trajectory  that  switches  from 
equation  (5)  to  (6)  or  vice  versa  at  9  =  45  degrees. 

From  observing  Figures  3  and  4,  the  sets  of  initial  conditions  yielding  complete  ECAV  trajectories  for  a 
full  90  degrees  (valid  initial  conditions)  are  small  compared  to  the  set  [0, 0.707],  which  contains  all  the  initial 
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Using  Eq.  (5) 


Using  Eq.  (5)  and  then  (6) 


-  Phantom  track,  a  =  2 

ECAV  trajectory 
_ r  0_,  =  1  boundary 

dot  1 

-  Radar  LOS 


Figure  3.  Constant-speed  ECAV  trajectory  solutions  for  a  constant-speed,  constant-course  phantom  track, 
starting  with  equation  (5). 


Using  Eq.  (6) 


Using  Eq.  (6)  and  then  (5) 


-  Phantom  track,  a  =  2 

ECAV  trajectory 
_ r-0 =  1  boundary 

dot  1 

-  Radar  LOS 


Figure  4.  Constant-speed  ECAV  trajectory  solutions  for  a  constant-speed,  constant-course  phantom  track, 
starting  with  equation  (6). 
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conditions  an  ECAV  could  start  deceiving  the  radar  from  for  a  =  2.  In  addition,  these  results  are  for  no 
additional  limitations  on  the  ECAV  dynamics.  The  set  size  of  these  valid  initial  conditions  varies  predictably 
with  both  phantom  track  speed  and  heading,  but  this  observation  is  not  pursued.  The  main  point  here  is 
that  the  constant-speed  limitations  on  the  ECAV  and  phantom  track  must  be  relaxed  to  a  realistic  range  to 
allow  the  ECAV  more  flexibility  in  choosing  its  starting  point  and  trajectory  in  the  deception  process. 


IV.  General  Theory  for  ECAV  Trajectory  Bounds  and  Solutions 

Because  the  assumption  of  a  constant-speed  ECAV  and  phantom  track  puts  severe  limitations  on  the 
ECAV  valid  initial  conditions  and  trajectories,  it  is  beneficial  to  move  the  scenario  closer  to  reality  by 
allowing  the  ECAV  and  phantom  track  speeds  to  vary  within  some  bounded  range.  In  addition,  it  is  now 
assumed  that  the  ECAV  has  two  fixed  antennas — each  mounted  on  one  side  of  the  ECAV  with  less  than  a 
90-degree  look-angle  from  its  central  axis.  Hence,  the  ECAV  can  no  longer  send  pulse  returns  to  a  radar  at 
any  course  angle,  <f>E- 

Mainlobe  deception  is  again  assumed,  so  the  ECAV  has  one  DOF,  which  will  be  constrained  in  dif¬ 
ferent  ways  to  create  several  ECAV  dynamic  systems  for  exploration  of  the  ECAV’s  flyable  range  given  a 
phantom  track.  For  ease  of  comparison,  the  following  non-dimensional  variables  are  used,  with  i>Em  =  1 
(average/nominal  ECAV  speed),  R0  =  1,  and  90  =  0. 


ve 


VEm  , 

Vt 
V Em 

VE 


VEn 


R 


Ro 

R 

Ro 


For  the  rest  of  this  study,  the  following  physically  realistic  ranges  are  used  in  determining  ECAV  trajectory 
bounds  and  solutions. 


ve  r 

-  ±20% 

(f)E  - 

-  90°  ±  60' 

a  r- 

-  ±20% 

Given  the  definitions  and  the  ranges  for  i>e,  <t> e >  and  ol  listed  above,  the  minimum  and  maximum  values  for 
these  same  variables  are  as  follows. 


VEmin 

=  0.8 

4*  Emin 

^  l<o 

II 

Vt-min  0.8a 

VEmax 

=  1.2 

(pEmax 

II 

05  Sf 

Vimax  ~  l-2a 

(12) 


Now  for  a  generic  phantom  track  and  the  previously  defined  non-dimensional  variables,  the  appropriate 
variables  and  their  relations  for  the  track  and  for  developing  ECAV  dynamic  systems  are  shown  below  in 
figure  5.  Definitions  of  the  basic  ECAV  and  phantom  track  variables  are  also  given  on  the  right  side  of  this 
figure. 

The  basic  equations  of  motion  for  the  ECAV  and  phantom  target,  using  polar  coordinates,  are  as  follows. 
Every  derivative  is  taken  with  respect  to  time.  An  equation  for  ECAV  time  delay,  td,  is  also  included. 


r  =  vE  cos  <))e ,  r(0)  =  r0 

(13) 

9=  —  sine/)#,  0(0)  =  0 
r 

(14) 

R=acos(/>T,  R(0)  =  1 

(15) 

0  =  —  sin  cf>T ,  0(0)  =  0 

r 

(16) 

td(t)  =  ^(R(t)-r(t)) 

(17) 
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a 


E  ECAV 

ve  ECAV  speed 

4>e  ECAV  course 

r  ECAV  range  from  radar 

9  ECAV/phantom  target 
bearing  from  radar 

T  phantom  target 

R  phantom  target  range 
from  radar 

vt  phantom  target  speed 


Figure  5.  Variables  and  relations  for  an  ECAV  with  variable  speed  and  a  0T  phantom  target  course 
generic  phantom  track  with  variable  speed. 


Without  loss  of  generality,  9  is  assumed  positive  and  given  by  some  pre-determined  phantom  track.  Six 
ECAV  dynamic  systems  are  now  presented  below,  which  can  each  be  developed  from  equations  (13)  and  (14) 
by  setting  one  ECAV  variable.  All  systems  include  an  algorithm  for  choosing  initial  conditions  that  satisfy 
the  ranges  in  (12),  ODEs  (ordinary  differential  equations)  necessary  to  solve  for  the  ECAV  trajectory,  and 
bounds  resulting  from  (12)  that  must  be  satisfied  for  all  time.  A  subscript  c  indicates  that  the  corresponding 
variable  is  constant.  Remember  that  9  is  a  function  of  time. 

A.  Constant  Speed  System  (vec) 

For  the  initial  conditions,  choose  Vec  and  r o  in  sequence  to  satisfy  the  following  initial  constraints. 


'UEmin  Iff  V Ec  if:  ^ Emax 


vEc  sin  cpE 


mm  ^  ^  ^Ec 

-  <r  o  < 

90  90 


Then  the  differential  equation  with  (+)  and/or  (— )  to  be  solved  is  given  as 

+  \lvEc2-{rO)2 

v -  r(0)  =  r0. 


- yvEc 2  -  M)2 

where  for  the  solution  r(f),  the  course  constraint 

r[t)9 


sin  (j>Emin  ^ 


VEc 


Vt 


(18) 


must  be  satisfied,  and  the  constraint 

s(t)  :=  i ’ec2  -  ( r(t)9 )2  >0  Vt 

must  also  be  satisfied;  however,  if  s(t)  =  0  and  condition  (8)  is  satisfied,  then  switching  to  solving  the  other 
half  of  equation  (18)  is  valid. 
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B.  Constant  Course  System  (4>ec) 

For  the  initial  conditions,  choose  <f>Ec  and  tq  in  sequence  to  satisfy  the  following  initial  constraints. 

pEmin  —  4* Ec  5s  pEmax 
'UEmin  Sill  (pEc  _  ^ Emax  Sin  (pEc 

—  <  r0  <  - 


9o  Sq 

Then  the  differential  equation  to  be  solved  is  given  as 

r  =  r9  cot  4>ec,  r(0)  =  r0, 

where  for  the  solution  r(t),  the  speed  constraint 

r(t)9 

VEmin  _  .  ,  ^  V Emax  vt 

Sin  CpEc 

must  be  satisfied. 

C.  Constant  Heading  System  (JiEc) 

For  the  initial  conditions,  choose  Hec  and  ro  in  sequence  to  satisfy  the  following  initial  constraints. 

< pEmin  +  9q  <  h-Ec  <  <pE  max  H"~  #0 
VEmin  Sin  (flEc  ^o)  _  V Emax  Sin  {h>Ec  ^o) 

Oo  do 

Then  the  differential  equation  to  be  solved  is  given  as 

f  =  r9cot(hEc  -  9(t)),  r(0)=ro, 

where  for  the  solution  r(t),  the  speed  constraint 

r(t)9 


(19) 


(20) 


V Emin  ^ 


sin  <pE{t) 


PP  VEmax  Vt 


must  be  satisfied. 


D.  Constant  Speed  Rate  System  (vec) 

For  the  initial  conditions,  choose  i)Ec,  Veo,  and  ro  in  sequence  to  satisfy  the  following  initial  constraints. 

VEmin  P4:  ^ Ec  ^  'S Emax 

VEmin  HP  V E0  ^  ^ Emax 
Veo  Sin  (p Emin  ,  ,  VEO 

- ■ -  <  r0  < 

do  90 

Then  the  differential  equation  to  be  solved  is  given  as 

—r9(r9  +  r9)  +  vE(t)vEc  r(0)  =  ro 


(0)  =  f0  =  ^ JvE0 2  -  {r09o)2 


(21) 


'  r(U)=rn  = 

where  for  the  solution  r(f),  the  speed  constraint 

V Emin  ~  ^e{4)  HP  V Emax  ^ 

must  be  satisfied,  the  course  constraint 

r(t)9 

Sin  (pEmin  —  7T7  V  t 

VE{t) 

must  be  satisfied,  and  the  constraint 

s(t)  :=  |r(i)|  >0  Vt 

must  also  be  satisfied;  however,  if  s(t )  =  0,  then  switching  to  — vEc  and  perturbing  f  across  zero  allows  the 
solution  of  equation  (21)  to  continue. 
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E.  Constant  Turn  Rate  System  ( Hec ) 

For  the  initial  conditions,  choose  h,Ec  and  then  (f>Eo  and  ro  in  sequence  to  satisfy  the  following  initial 
constraints. 

$Emin  ^  $E0  ^Emax 
V Emin  sin  0  EO  _  V Emax  sin  0  EO 

—  <r0<  - 


On 


On 


Then  the  differential  equation  to  be  solved  is  given  as 

„  _  f(rO  +  rO)  -  VE2{t)(h.Ec  -  0)  r(0)  =  r0 

rO  ’  r(0)  =f0  =  ro0o  cot  fy eo 

where  for  the  solution  r(t),  the  speed  constraint 

VEmin  —  Ve  (f  )  ^  V Emax 

must  be  satisfied,  and  the  course  constraint 

r{t)6 

sin  (p Emin  —  77T  Vt 


(22) 


Ve(*) 


must  be  satisfied. 


F.  Constant  Acceleration  System  (aEc) 

For  the  initial  conditions,  choose  agc  and  then  veo  and  ro  in  sequence  to  satisfy  the  following  initial  con¬ 
straints. 

VEmin  VEO  —  VEmax 


veo  sin  <f>E  r i 

Jo 


,  ,  VEO 

-ro-17 


Then  the  differential  equation  with  (+)  and/or  (— )  parts  to  be  solved  is  given  as 


r 02  +  y  ciEc2  —  {rO  +  2  rO)2 
r6 2  —  aEc2  ^  (vO  +  2 rO)2 


r( 0)  =  r0 


f(0)  =  r0  =  Y  veo2  -  {roOo)2 

where  for  the  solution  r(f),  the  speed  constraint 

VEmin  —  Ve{1)  ^  VEmax  Vt 

must  be  satisfied,  the  course  constraint 

r(t)6 

Sin  (pEmin  —  7T7  V  t 

vsyt) 

must  be  satisfied,  and  the  constraint 

s(t)  :=  aEc 2  —  {r(t)9  +  2 r0)2  >0  Vf 


(23) 


must  also  be  satisfied. 

System  A  above  simply  contains  a  reformulation  of  equations  (5)  and  (6)  for  a  range  of  constant  ECAV 
speeds.  System  B  allows  the  choice  of  a  constant  ECAV  course,  which  is  actually  relative  to  6  (see  figure  5 
above);  a  course  of  90  degrees  would  result  in  a  circular  trajectory.  Equation  (19)  in  system  B  may  actually 
be  solved  explicitly  for  r(f),  but  this  solution  is  not  a  concern.  System  C  allows  the  choice  of  a  constant  ECAV 
heading,  which  means  that  the  ECAV  will  fly  in  a  straight  line  with  the  specified  heading.  Equation  (20) 
may  also  be  solved  explicitly  for  r(t)  if  the  system  is  autonomous,  i.e.  if  6  is  not  a  function  of  time.  System  D 
allows  the  choice  of  a  constant  ECAV  speed  rate  and  is  not  too  useful  since  s(t)  quickly  reaches  zero  for 
most  initial  conditions.  System  E  allows  the  choice  of  a  constant  ECAV  turn  rate  and  is  extremely  useful 
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because  it  can  often  be  used  to  find  the  minimum  and  maximum  initial  conditions  for  which  a  flyable  ECAV 
trajectory  exists  for  a  specified  range  of  9.  System  F  allows  the  choice  of  a  constant  ECAV  acceleration  (time 
derivative  of  ECAV  velocity)  and  is  more  useful  when  a  given  phantom  track  is  circular  in  form  so  that  a 
constant  acceleration  greater  (less)  than  the  inherent  centripetal  acceleration  causes  the  ECAV  to  spiral  in 
(out)  relative  to  the  radar. 

The  variable  definitions  and  minimum/maximum  values  for  ECAV  speed  and  course  in  (12)  give  rise 
to  the  following  r0-boundary  lines  for  the  ECAV,  which  are  similar  in  concept  to  condition  (7)  with  a 
constant-speed  ECAV. 


r(t)  = 

V Emin  0*8 

(24) 

9  9 

r(t)  = 

VEmax  1*2 

(25) 

9  9 

r(t)  = 

V Emin  Sin  4* Emin 

0.4 

(26) 

9 

9 

r(t)  = 

V Emax  Sin  (ft  Emin 

0.6 

(27) 

9 

9 

The  ECAV’s  range  may  not  exceed  boundary  (25)  at  any  time  because  this  would  require  it  to  fly  faster 
than  its  maximum  speed.  If  the  ECAV  range  decreases  below  boundary  (24),  the  ECAV  may  continue  its 
trajectory,  but  may  not  fly  perpendicular  to  its  current  LOS  since  it  must  have  a  nonzero  r-component  to  fly 
above  its  minimum  speed.  The  ECAV’s  range  may  not  decrease  below  boundary  (26)  at  any  time  because 
this  would  rotate  its  fixed  antennas  out  of  range  of  the  radar  even  at  minimum  speed.  If  the  ECAV  range 
decreases  below  boundary  (27),  the  ECAV  may  continue  its  trajectory,  but  only  at  a  speed  less  than  its 
maximum  speed. 

For  this  analysis  on  ECAV  trajectory  bounds,  it  is  beneficial  to  convert  the  phantom  target  speed  range 
into  a  larger  pseudo-range  for  the  ECAV  speed  with  the  phantom  target  speed  resumed  to  a  virtual  constant. 
To  make  this  conversion,  the  effect  on  the  r0-boundary  lines  (24)-(27)  due  to  varying  a  down  to  amin  and  up 
to  amax  is  determined  and  maximized  to  produce  a  new  set  of  r0-boundary  lines  for  the  ECAV  speed  pseudo¬ 
range.  Hence,  the  effect  of  a  on  6  must  first  be  determined  for  various  phantom  tracks.  The  dependence  of 
9  on  a  for  a  straight  phantom  track  is  given  as  follows,  using  equations  (9)— (11) ,  with  9  solved  explicitly  as 
a  function  of  9. 

m  =  — - ^ -  (28) 

1+  (  .  ,sinfl  )  —  2 cos if  .  ,sine 

ysin(7 T—'ip—O)  I  •  sin(7r— yj— 6) 


Conveniently,  9  is  directly  proportional  to  a  for  any  given  value  of  9 ,  and  this  result  can  graphically  be 
shown  to  hold  true  also  when  the  phantom  track  is  circular  and  referenced  to  a  radar  not  at  its  center  as 
in  equations  (42)-(44)  below.  Therefore,  the  new  rfl-boundary  lines  are  as  follows  for  a  straight  or  circular 
phantom  track,  where  VEpsmin  and  VEpsmax  are  defined  in  equations  (29)  and  (30),  respectively. 


r{t)  = 

VEmin 

VEpsmin 

9 

ipmax  Ot)0 

r(t )  = 

VEmax 

VEpsmax 

9 

(q(  OLmin)@ 

r(t)  = 

VEpsmin  S1H  (ft Emin 

0.33 

9 

9 

r(f)  = 

VEpsmax  Sin  (ft Emin 

0.75 

9 

9 

0.67 

9 

(29) 

1.5 

(30) 

T 

(31) 

(32) 

Figure  6  below  is  an  exported  frame  from  a  MATLAB-generated  video  and  illustrates  visually  how  rO- 
boundaries  (24)  and  (25)  are  changed  to  r0-boundaries  (29)  and  (30),  respectively,  by  allowing  varia¬ 
tion  in  a  to  move  the  ECAV  speed  range  in  or  out  to  further  limits.  This  visual  explanation  is  simi¬ 
lar  for  the  changes  made  to  convert  boundaries  (26)  and  (27)  to  boundaries  (31)  and  (32),  respectively. 
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Phantom  track,  a  =  2  +  0% 


r-Bd  boundary  for  v£  =  1 
r-Qd  boundary  range  for  x>F±  20% 


r 0  .  boundary  for 

dot  Epsmax 


E~ 

=  1. 

=  O.f 


In  figure  6,  the  gray  boundary  represents  the 
outer  limit  where  a  is  decreased  by  20%  to  amin, 
which  effectively  moves  the  black  boundary  range 
for  ve  ±  20%  (the  area  between  the  black  curves) 
out  to  the  gray  boundary;  at  this  point,  the  outer 
black  boundary  coincides  with  the  gray  boundary. 

Likewise,  the  cyan  boundary  represents  the  inner 
limit  where  the  black  boundary  range  may  be  moved 
when  a  is  increased  by  20%  to  amax.  It  is  impor¬ 
tant  to  realize  that  these  ultimate  lower  and  upper 
r(9-boundaries  (the  cyan  and  gray  curves)  may  not 
coexist  at  any  instant  in  time  because  a  may  only 
be  one  value  at  any  time.  This  fact  is  important  pri¬ 
marily  in  the  case  of  multiple  ECAVs  and  multiple 
radars  because  it  constrains  all  ECAVs  to  be  within 
the  ^-boundary  range  (the  area  between  the  black 
curves)  at  any  given  time,  even  though  that  range 
may  move  in  or  out  as  shown  in  figure  6.  When 
changing  a  to  accommodate  an  ECAV’s  desired  po¬ 
sition  or  speed,  a  will  have  to  modulate  in  a  contin¬ 
uous  and  believable  fashion,  especially  if  the  radar 
is  operating  at  a  high  PRF. 

To  further  enhance  the  assessment  of  an  ECAV’s  flyable  range  given  a  phantom  track  and  some  range  of 
9 ,  the  following  speed  isolines  may  also  be  used,  where  R(t)  is  specified  by  the  phantom  track  being  created. 


Black  boundary 
curves  move 
•  together  in 
one  direction 


Radar 

Figure  6.  Conversion  of  phantom  track  speed  range 
into  larger  ECAV  speed  pseudo-range. 


r(f)  =  VErsminR(t) 
a 

(33) 

r(t)  -  VEpamaxR(t) 

(34) 

V.  ECAV  Bounds  for  a  Straight  Phantom  Track 

THE  theory  presented  in  the  previous  section  for  treating  ranges  on  ECAV  speed,  course  (region  of  antenna 
operation) ,  and  the  speed  of  the  phantom  track  is  now  utilized  to  conduct  a  survey  of  the  ECAV  bounds 
for  a  constant-course  phantom  track  through  the  solution  of  ordinary  differential  equations.  The  parametric 
equations  for  this  straight  phantom  track,  taken  from  Ref.  2,  are  shown  below  in  polar  coordinates  and  are 
identical  to  equations  (9)-(ll). 


R(t)  =  \/Y+~c^t^~^2at-  cos  ip 
'  at  sin  ip N 

W 

a  sin  'll) 


6{t)  =  arcsin 

m  = 


1  +  a2t2  —  2 at  cos  ij) 


(35) 

(36) 

(37) 


Using  equations  (35)-(37)  and  with  a  =  2  and  if)  =  45  degrees,  a  phantom  track  is  plotted  along  with 
the  boundaries  and  isolines  corresponding  to  (29)-(32)  and  (33)-(34),  respectively,  in  figure  7  below.  The 
yellow  lines  represent  a  flyable  ECAV  range  for  9  =  90  degrees.  This  flyable  range  represents  the  union  of  all 
positions  the  ECAV  could  visit  on  certain  trajectories  and  still  be  able  to  create  a  phantom  track  through 
90  degrees. 

Many  ECAV  trajectories — solved  using  the  six  dynamic  systems  A  through  F  presented  in  the  previous 
section — were  used  to  test  and  verify  the  results  shown  in  figure  7.  The  flyable  range  bounded  by  the 
pseudo-speed  isolines  is  valid  since  any  trajectory  parallel  to  and  between  these  two  isolines  is  within  the 
ECAV  speed  pseudo-range  and  course  range.  The  additional  part  of  the  flyable  range,  defined  where  the 
green  curves  are  less  than  the  ve  =  0.67  isoline,  is  actually  bounded  by  a  (minimum)  constant-speed  ECAV 
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Phantom  track,  a  =  2 
v  =  1.5  isoline 

v£  =  0.67  isoline 
r-0dof=  1.5  boundary 
r-Qdo[  =  0.67  boundary 
e  =  30P( ue=1.5) 

^  =  30°  (ve=0.67) 

Radar  LOS 

Bound  defined  by  ECAV  traj. 


0  =  90 


\|/  =  45° 


Figure  7.  ECAV  flyable  range  for  a  constant-course  phantom  track,  6  =  90  degrees. 


trajectory  solution  using  system  A,  where  the  ECAV  is  initially /finally  on  the  red  boundary,  i.e.  its  course 
is  at  4>Emin  for  both  0  =  0  and  0  =  90  degrees.  For  an  ECAV  starting  on  the  right  trajectory  curve  and 
running  into  the  r9  =  0.67  boundary,  it  could  then  switch  to  a  different  system — such  as  system  B  with 
(f>Ec  =  90  degrees — to  create  a  circular  trajectory  and  then  switch  back  to  the  constant  speed  system  with 
vec  =  0.67  when  it  again  reached  the  boundary.  This  line  of  reasoning  as  well  as  many  flyable  ECAV 
trajectories  produced  using  a  constant  turn  rate  with  system  E  provide  ample  evidence  for  the  flyable  range 
below  the  speed  isolines  in  figure  7. 

The  right  “black  hole”  in  figure  7  is  an  area  that  can  be  entered  by  the  ECAV,  but  once  the  ECAV  is  in 
this  region,  it  cannot  exit  back  into  the  flyable  range;  it  will  run  into  the  rO  =  1.5  boundary  and  stop.  To 
explain  this  mathematically,  the  following  equation  for  the  ECAV  speed  at  any  time  instance  is  used. 


ve  =  \J  r2  +  (r0)2  (38) 

Taking  a  point  on  the  Ve  =  1-5  isoline  that  bounds  the  black  hole,  if  the  absolute  value  of  r  is  decreased, 
which  corresponds  to  entering  the  black  hole,  then  equation  (38)  says  that  ve  will  also  decrease.  This  allows 
the  ECAV  to  enter  the  black  hole  because  in  doing  so  it  will  be  flying  within  its  speed  pseudo-range  or  below 
VEpsmax-  However,  for  an  ECAV  to  exit,  the  absolute  value  of  f  would  have  to  increase,  which  by  equation 
(38)  requires  that  the  ECAV  fly  faster  than  its  maximum  pseudo-speed  of  1.5.  A  similar  argument  can  be 
used  to  show  that  the  region  labeled  “off  limits”  is  an  area  that  can  be  exited  by  the  ECAV  but  not  entered 
from  its  flyable  range.  The  left  “black  hole”  is  an  area  that,  once  entered,  will  require  that  the  ECAV  turn  in 
towards  the  radar  at  a  progressive  rate  to  stay  within  its  speed  pseudo-range;  at  some  point  in  this  process, 
the  ECAV  will  reach  its  course  bound  between  the  orange  and  red  curves  where  its  fixed  antennas  rotate 
out  of  range  of  the  radar — all  before  0  =  90  degrees  is  reached. 

In  addition  to  the  speed  and  course  bounds  discussed  above,  there  exist  ultimate  0-bounds  on  ECAV 
trajectories  given  a  straight  phantom  track.  For  an  ECAV  flying  parallel  to  a  speed  isoline,  it  would  reach 
its  course  bound  or  antenna  limit  at  60  degrees  off  of  the  shortest  LOS  from  the  radar  to  the  phantom  track. 
This  corresponds  directly  to  the  range  chosen  for  (f>E  in  the  previous  section,  which  is  90  ±  60  degrees.  It 
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may  be  possible  for  the  ECAV  to  continue  the  phantom  track  for  slightly  greater  values  of  9  by  starting  on 
a  high  speed  isoline  and  turning  in  towards  the  radar  right  before  it  reaches  its  course  bound;  however,  this 
possibility  was  not  explored. 

The  analysis  in  this  section  is  modular  in  that  it  is  applicable  to  any  number  of  ECAVs  and  the  same 
number  of  radars,  each  on  an  individual  basis.  The  parameter  ij)  of  the  phantom  track  is  simply  calculated 
for  each  radar,  and  a  different  set  of  modified  phantom  track  equations  are  then  used  for  each  radar/ECAV 
pair  so  that  together  the  team  creates  one  coherent  straight  phantom  track  with  speed  a. 


VI.  ECAV  Bounds  for  a  Simple  Circular  Phantom  Track 

The  general  theory  for  ECAV  trajectory  bounds  and  solutions  is  again  utilized  to  conduct  a  survey  of  the 
ECAV  bounds  for  a  circular  phantom  track  with  the  radar  placed  at  the  circle’s  center.  The  equations 
for  this  type  of  phantom  track,  taken  from  Ref.  2,  are  shown  below  in  polar  coordinates. 


R(t)  =  1 

(39) 

9(t)  =  at 

(40) 

9{i)  =  a 

(41) 

With  equations  (39)-(41)  representing  the  circular  phantom  track,  ECAV  systems  A  through  F  are  now 
autonomous  and  may  be  analyzed  in  the  phase  plane  if  so  desired.  Because  of  this  simplification,  the 
simple  circular  phantom  track  presents  a  good  opportunity  to  more  thoroughly  analyze  the  ECAV  dynamic 
systems  A  through  F.  Using  equations  (39)-(41)  and  with  a  =  2,  a  phantom  track  is  plotted  along  with  the 
boundaries  corresponding  to  (29)-(32)  in  figure  8  below.  The  yellow  lines  represent  a  flyable  ECAV  range 
for  6  =  180  degrees.  This  flyable  range  represents  the  union  of  all  positions  the  ECAV  could  visit  on  certain 
trajectories  and  still  be  able  to  create  a  phantom  track  through  180  degrees. 


-  Phantom  track,  a  =  2 

r-Qdot  =  1.5  boundary 

re  =  0.67  boundary 
§E  =  30°  (ue=1.5) 

-  ty£=  30°  (v£=  0.67) 

-  Radar  LOS 

-  Bound  defined  by  ECAV  traj. 


Figure  8.  ECAV  flyable  range  for  a  simple  circular  phantom  track  and  9  =  180  degrees. 

Many  ECAV  trajectories — solved  using  the  six  dynamic  systems  A  through  F  presented  earlier — were 
used  to  test  and  verify  the  results  shown  in  figure  8.  For  the  simple  circular  phantom  track,  the  minimum 
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and  maximum  pseudo-speed  isolines  are  not  identified  because  they  coincide  with  the  rQ  =  0.67  and  r9  =  1.5 
boundaries,  respectively;  this  is  advantageous  because  it  means  that  the  ECAV  has  no  ultimate  0-bounds 
as  with  the  straight  phantom  track.  Justification  for  the  flyable  range  presented  above  in  figure  8  is  not 
presented  here  because  the  simple  circular  phantom  track  is  actually  a  special  case  of  the  general  circular 
phantom  track,  which  is  presented  in  the  next  section. 

The  region  labeled  “off  limits”  in  figure  8  is  an  area  that  cannot  be  entered  or  exited  by  the  ECAV  because 
it  is  bounded  by  r9  =  1.5.  Taking  a  point  on  this  boundary,  to  enter  or  exit  the  off  limits  region  requires 
r  to  be  nonzero,  which  increases  the  ECAV  speed  above  its  pseudo-maximum  according  to  equation  (38). 
The  “black  hole”  is  an  area  that,  once  entered,  will  require  that  the  ECAV  turn  in  towards  the  radar  at  a 
progressive  rate  to  stay  within  its  speed  pseudo-range;  at  some  point  in  this  process,  the  ECAV  will  reach 
its  course  bound  between  the  orange  and  red  curves  where  its  fixed  antennas  rotate  out  of  range  of  the 
radar — all  before  0  =  180  degrees  is  reached. 

Although  simple  and  useful  for  understanding  ECAV  systems  and  trajectory  bounds,  the  analysis  in  this 
section  alone  is  not  modular.  It  is  applicable  to  only  one  ECAV  trying  to  deceive  a  radar  located  at  the 
center  of  the  circular  phantom  track.  However,  this  information  combined  with  the  general  circular  phantom 
track  analysis  presented  in  the  next  section  will  provide  the  modularity  needed  for  a  team  of  ECAVs  to 
create  one  coherent  circular  phantom  track  for  a  radar  network. 


VII.  ECAV  Bounds  for  a  General  Circular  Phantom  Track 


The  general  theory  for  ECAV  trajectory  bounds  and  solutions  is  once  more  utilized  to  conduct  a  survey 
of  the  ECAV  bounds  for  a  circular  phantom  track  with  the  radar  placed  arbitrarily  with  reference  to 
the  circle’s  center.  The  equations  for  this  type  of  phantom  track  are  shown  below  in  polar  coordinates.  Only 
two  of  the  phantom  track’s  three  additional  parameters,  a,  6,  and  p,  must  be  specified;  a  and  b  represent  the 
horizontal  and  vertical  distances  from  the  circle  center  to  the  radar. 


I  O'  CX 

R(t)  =  ,  1  +  2  ap  —  2apcos  —  t  —  26psin  —  t 
V  P  P 

6{t)  =  sgn  ^  (x(t)  -7 r)  +  £,  7r-£<f<7r  +  £ 

m=  a 


,  a  a 

i  p  -  acos  -t  -  osin  —t 

R2{t)  \  p  p 


(42) 

(43) 

(44) 


The  mathematical  definitions  for  \(t)  and  the  parameters  p  and  £  are  as  follows. 


/\  /  .  lx  7  u:  272  CX  .  CX 

X\t)  =  arg  \  p  a  sm  t  —  b  cos  —t,a  +  6"  —  ap  cos  —  t  —  bp  sin  —  t 

P  P 


p  =  a+  \J\  —  b2 
£  =  arg(6,  a) 

Using  equations  (42)-(44)  and  with  a  =  2,  a  =  —0.38,  and  b  =  0.22  (radar  placed  inside  the  circle),  a 
phantom  track  is  plotted  along  with  the  boundaries  and  isolines  corresponding  to  (29)-(32)  and  (33)-(34), 
respectively,  in  figure  9  below.  The  yellow  lines  represent  a  flyable  ECAV  range  for  a  circular  phantom  track 
of  180  degrees,  which  corresponds  to  9  =  227  degrees.  This  flyable  range  represents  the  union  of  all  positions 
the  ECAV  could  visit  on  certain  trajectories  and  still  be  able  to  create  a  180-degree  circular  phantom  track. 

Many  ECAV  trajectories — solved  using  dynamic  systems  A  and  E  presented  earlier — were  used  to  test 
and  verify  the  results  shown  in  figure  9.  For  the  general  circular  phantom  track,  the  relevant  isolines  and 
boundaries  plotted  in  figure  9  look  similar  to  those  plotted  for  the  straight  phantom  track  (see  figure  7  above) 
except  curved  around  in  an  egg-like  shape.  The  flyable  range  bounded  by  the  speed  isolines  is  valid  since  any 
ECAV  trajectory  coinciding  with  a  u^-isoline  between  these  two  will  fly  at  a  constant  speed  and  be  within 
the  ECAV  speed  pseudo-range  and  course  range.  The  additional  part  of  the  flyable  range,  defined  where 
the  green  curves  are  less  than  the  i>e  =  0.67  isoline,  is  actually  bounded  by  a  (minimum)  constant-speed 
ECAV  trajectory  solution  using  system  A,  where  the  ECAV  is  initially /finally  on  the  red  boundary,  i.e.  its 
course  is  at  (j>Emin  for  both  0  =  0  and  0  =  90  degrees.  For  an  ECAV  starting  on  the  right  trajectory  curve 
and  running  into  the  r0  =  0.67  boundary,  it  could  then  turn  into  an  appropriate  speed  isoline  in  between 
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-  Phantom  track,  a  =  2 

\>  =  1.5  isoline 

\>E  =  0.67  isoline 
r'®dot  =  1  5  boundary 
r-0  =  0.67  boundary 

=  30°  (u£  =  1 .5) 

■-  =  30°  (u£  =  0.67) 

-  Radar  LOS 

-  Bound  defined  by  ECAV  traj. 


Figure  9.  ECAV  flyable  range  for  a  general  180-degree  circular  phantom  track  with  the  radar  placed  inside 
the  circle. 


the  minimum  and  maximum  isolines  and  switch  back  to  the  constant  velocity  system  with  vec  —  0.67  when 
it  again  reached  the  boundary.  This  line  of  reasoning  as  well  as  some  flyable  ECAV  trajectories  produced 
using  a  constant  turn  rate  with  system  E  provide  evidence  for  the  flyable  range  below  the  speed  isolines  in 
figure  9.  It  has  been  verified  that  mirroring  the  radar’s  location  about  the  phantom  circle  center’s  y-axis 
(vertical  axis)  produces  symmetric  bounds  and  flyable  ECAV  trajectories  about  this  axis. 

The  “black  hole”  and  “off  limits”  regions  are  analogous  to  those  shown  for  the  straight  phantom  track 
in  figure  7  above.  Their  explanations  are  likewise  similar  and  so  are  omitted. 

Note  that  choosing  the  radar  outside  the  circular  path  of  the  phantom  track  will  severely  limit  how  long 
in  time  or  9  the  phantom  track  may  be  generated.  When  the  phantom  track  heading  comes  close  to  pointing 
towards  the  radar  ( 6  close  to  zero),  creation  of  the  track  will  stop.  This  result  is  due  to  the  course/antenna 
bounds  imposed  on  the  ECAV. 

The  analysis  in  this  section,  combined  with  that  for  the  simple  circular  phantom  track  in  the  previous 
section,  now  provides  a  modular  approach  for  a  circular  phantom  track  applicable  to  any  number  of  ECAVs 
and  the  same  number  of  radars,  each  on  an  individual  basis.  The  parameters  a  and  b  of  the  phantom  track 
are  simply  recalculated  for  each  radar,  and  a  different  set  of  modified  phantom  track  equations  are  then  used 
for  each  radar/EC AV  pair  so  that  together  the  team  creates  one  coherent  circular  phantom  track  with  speed 
a.  However,  if  any  radar  is  placed  outside  the  circular  phantom  track,  the  track  will  be  generated  for  less 
than  180  degrees  due  to  the  ECAVs’  course/antenna  limitations  from  (12). 
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VIII.  General  Bounds  for  ECAV  Initial  Conditions  and  Flyable  Ranges 


Based  on  the  ECAV  flyable  ranges  developed  for  straight  and  circular  phantom  tracks  and  presented 
in  figures  7  through  9  above,  the  following  equations  define  and  summarize  the  range  of  valid  initial 
conditions  for  an  individual  ith  ECAV  to  start  from  and  generate  a  straight  or  circular  phantom  track  for  a 
desired  range  of  9 ,  where  i?0i  is  a  dimensional  length. 


(fOmin)i 

VEpsmin  Sin  4* Emin  j-, 

—  ■  -Ko  i 

(45) 

9oi 

{j’O  max)i 

( V Epsmax  \  t-> 

=  mm  ^  ,  1J  Roi 

(46) 

Due  to  the  conversion  of  the  phantom  track  speed  range  into  a  larger  ECAV  speed  pseudo-range,  the  initial 
conditions  of  multiple  ECAVs  are  constrained  to  stay  within  the  original  r0-boundary  range  corresponding 
to  their  specified  UE-range  of  ±20%  for  whatever  value  a  actually  is  at  any  time  instance  (see  figure  6  above). 
Therefore,  equations  (45)  and  (46)  specify  a  maximum  range,  in  which  the  selected  initial  condition  values 
for  each  ECAV  must  also  satisfy  the  following  conditions,  where  I  contains  an  ordered  set  of  index  numbers 
for  each  ECAV  participating. 


if 

then 

if 

then 


.  'UEmin  , , 
TOi  ±  Roi 


r0i  <  mm 
jei 


r0j 


Rt 


■0  j 


.  VEmax  T) 
T'Oi  ^  -*M) i 


a 

r^i  >  max 
jei 


(  r0j 

\Roj 


for  any  i  £  I 

VEmax  D  w  •  _  r 

- R0i  vi  £  I 

V Emin 


for  any  i  £  I 


V Emin 
'VEmax 


-R'Oi 


Mi  £  I 


(47) 


(48) 


In  addition  to  the  range  of  valid  initial  conditions,  the  time-dependent  upper  and  lower  bounds  on  the 
flyable  range  of  an  individual  ith  ECAV  generating  a  straight  or  circular  phantom  track  for  some  range  of  9 
are  given  below,  where  Ri{t)  is  dimensional. 

(rmm),(t)  =  min  (  VEp^nin  R4  ( t ) ,  rf 0  (t) ,  rf f  (t)  j  (49) 

(■ rmax)i{t )  =  min  ( VEps™ax  ,  l)  Rt(t)  (50) 

In  equation  (49),  the  solutions  rf°(t)  and  rf*{t)  are  to  equation  (18)  from  ECAV  system  A,  and  both  use 
VEpsmin  for  their  constant  speed.  In  addition,  rf°(t)  uses  (romin)i  for  its  initial  value — determined  from 
equation  (45);  and  rf^{t)  uses  ( T’/mjn)i  f°r  its  final  value — determined  from  equation  (45)  with  subscript  / 
substituted  for  subscript  0.  These  special  bounding  solutions  are  represented  graphically  for  specific  cases  in 
figures  7  through  9  as  green  ECAV  trajectories.  Equations  (49)  and  (50)  specify  a  maximum  time-dependent 
flyable  range,  in  which  each  ECAV’s  position  must  also  satisfy  conditions  (47)  and  (48)  for  all  time,  with 
the  following  substitutions. 

n(t)  — >  r0i 
Ri(t)  — >  Roi 


IX.  Formulation  of  a  Decentralized  Cooperative  Problem 

AS  a  first  cut,  the  decentralized  cooperative  control  problem — one  of  team  coordination — is  posed  using 
the  concept  of  coordination  functions  and  the  general  bounds  in  (45)-(50);  the  aim  here  is  to  show  one 
way  the  results  from  the  previous  section  may  be  usefully  applied  to  control  the  ECAV  team.  This  problem 
is  illustrated  for  a  team  of  three  ECAVs  using  mainlobe  deception  on  three  radars  by  creating  a  straight 
phantom  track.  To  begin,  assume  a  constant-course  phantom  track  from  right  to  left  is  desired,  and  that 
bounded  ranges  for  ve,  &e,  and  a  have  already  been  given  as  in  (12).  Let  (3  be  the  absolute  heading  of  the 
phantom  track  and  n  be  the  number  of  radars — equal  to  the  number  of  ECAVs — then  mathematically  define 
the  following  parameters. 
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I  =  [l,n]  £  N 


i  =  1, 2, . . . ,  n 
rdi  =  ( rdxi,rdyi )  £  R2 
C'lTclCi^x'j  —  fdyi  d b  \J (Rmax')i‘^  (*£  f dxi 
COJZei  =(/3+^)±(^  -  4>Emin^j 


In  the  definitions  above,  the  location  of  each  radar  is  rdi.  The  circle  and  cone  equations  represent  subsets  of 
the  total  space  where  deception  is  possible  for  each  radar/ECAV  pair  given  the  radar’s  maximum  operational 
range  and  the  ECAV’s  course/antenna  limitations  from  (12).  With  the  above  parameter  definitions,  the 
following  function  can  be  defined,  which  maps  these  parameters  into  the  set  of  candidate  functions  that  will 
later  be  used  as  inputs  to  each  ECAV’s  coordination  function. 


hr(rdiei,  {Rmax)iel,  (3,  4>Emin)  *  Vi3 


The  set  consists  of  all  functions  f(x)  specifying  possible  phantom  tracks  that  all  ECAVs  can  deceptively 
create  for  the  participating  radar  network;  thus,  each  function  lies  in  the  space  B — common  to  all  circle 
and  cone  bounded  areas.  Mathematically,  each  function  f(x)  satisfies  the  following  conditions;  the  second 
condition  is  specific  to  the  straight  phantom  track  case. 


(*>/(*))  £  circlet 

) n  (q  conei ^  =:  B 


df  .  _ 

—  =  tan  fj 
ax 


Figure  10  below  illustrates  the  mathematics  introduced  so  far,  where  n  =  3  radars.  The  blue  lines 
represent  functions  in  and  are  all  possible  candidates  for  the  team-optimal  phantom  track,  i.e.  each  may 
be  generated  given  the  three-radar  situation  and  taking  into  account  each  radar’s  maximum  range  (gray 
solid  lines)  and  the  course/ antenna  bounds  associated  with  the  ECAVs  (black  dashed  lines). 

Three  variables  have  yet  to  be  chosen  before  a  given  phantom  track  is  decided  for  the  ECAV  team.  These 
variables  and  one  approach  to  choosing  them  are  given  below. 

f{x)  Let  the  phantom  track  as  a  function  of  x  be  determined  by  the  cooperative  control  process. 


Xq  For  simplicity,  assume  that  the  longest  length  of  phantom  track  is  desirable,  so  choose  Xq  by  using  the 
right-side  boundary  of  B,  Xq (y),  and  set  xo(y)  >  x  by  the  right-to-left  assumption. 


a  Since  the  phantom  target  speed  already  has  an  allowable  range  of  ±20%,  assume  its  nominal  or  mean 
value  is  fixed  by  the  type  of  phantom  target  being  created,  but  ensure  a  >  1.5  to  use  the  full  ECAV 
speed  and  phantom  target  speed  ranges. 

Given  the  variables  /( x),  Xq,  and  a,  the  following  function  may  be  defined,  which  uses  these  variables 
to  determine  valid  initial  condition  bounds  for  each  ECAV — using  equations  (45)  and  (46) — and/or  time- 
dependent  bounds  for  the  flyable  range  of  each  ECAV — using  equations  (49)  and  (50). 


hr(rdz,  f(x),x0,a) 


/ommdOmaxjj  ;  ^0 i 
min  (t)  5  ^ max  (*)]». 


The  intermediate  calculations  needed  for  using  equations  (45)  and  (46)  to  find  romin  and  romax  for  each 
ECAV  are  given  below. 


Roi  =  11(^0)  fixo))  -  rdi  ||2 

9oi  =  arg (f(x0)  -  rdyi ,  x0  -  rdxi) 

,  a  df 

i  —  arctan  — 
ax 

90i  =  a  sin  ipi 


18  of  20 


American  Institute  of  Aeronautics  and  Astronautics 


The  additional  intermediate  calculations  needed  for  using  equations  (49)  and  (50)  to  find  rmin(t)  and  rmax{t ) 
for  each  ECAV  are  given  below. 


Riit)  =  i?o?:\/l  +  o.2t 2  —  2  at  cos 
'  at  sin  ipj, 


Ri{t) 

a  sin  ipi 


6i{t )  =  arcsin 

Q.(f\  —  _ 

1  1  +  a2t2  —  2 at  cos 


Of  the  seven  intermediate  equations  listed  above,  the  last  five  are  specific  to  the  straight  phantom  track 
case.  Recall  that  conditions  (47)  and  (48)  must  be  applied  to  the  choice  of  initial  conditions  and  positions 
over  time  for  each  ECAV. 

Two  cost  functions  are  now  defined  for  each  ECAV  and  for  a  given  phantom  track  f(x)  as  follows,  where 
the  cost  could  be  figured  based  on  a  variety  of  factors  such  as  fuel  or  radiated  power  required  for  deception. 


icOStj  ( [/‘Omin  5  T'0max\ii  @0i) 
tcosti  ( [r  o  min  i  T'0max\ii  @0 it  \T  minify)  ^ max  (*)]<.*<(*)) 

The  initial  cost  of  each  ECAV,  icost ,  measures  how  hard  it  would  be  for  that  ECAV  to  move  from  its 
current  position  to  the  initial  position  required  for  commencing  a  given  phantom  track.  Note  that  the 
time-dependent  cost  of  each  ECAV,  tcost,  is  not  defined  to  depend  on  the  positions  of  other  ECAVs  despite 
conditions  (47)  and  (48);  rather,  the  cost  simply  depends  on  the  maximum  flyable  range  of  each  ECAV 
before  these  conditions  are  applied  to  actual  positions. 

Finally,  a  cost  based  on  the  initial  conditions  and  flyable  ranges  of  each  ECAV  can  be  computed  as 
a  function  of  the  track  choice,  f(x) — the  only  variable  left  to  be  chosen — which  leads  to  the  following 
coordination  function  for  each  ECAV  with  three  variations  possible. 

{icosti 
tcosti(t ) 

icosti  +  tcosti(t) 

Each  ECAV  passes  only  a  coordination  function,  hi,  to  the  team  leader.  The  team  leader  then  uses  all  n 
coordination  functions  to  find  the  function,  f(x),  that  will  minimize  the  total  cost  of  the  team  and  specify  the 
team-optimal  phantom  track  for  the  ECAVs  to  create.  If  the  first  version  of  the  coordination  function  is  used, 
then  each  ECAV  must  know  four  things:  its  own  radar’s  location  and  maximum  range,  its  initial  condition 
line  xq  (y),  and  the  phantom  track  speed.  If  the  second  or  third  versions  of  the  coordination  function  are  used, 
then  each  ECAV  must  additionally  know  the  locations  and  maximum  ranges  of  all  additional  radars  involved 
because  its  time-dependent  flyable  range  depends  on  where  the  phantom  track  stops  as  well  as  starts.  If  the 
flyable  ranges  were  defined  to  be  only  the  space  in  between  minimum  and  maximum  pseudo-speed  isolines 
(see  figures  7  through  9),  then  the  flyable  range  for  each  ECAV  would  only  be  a  function  of  x  and  not  t. 
Thus,  each  ECAV  would  only  need  to  know  its  own  radar’s  location  and  maximum  range  if  version  two  of 
the  coordination  function  above  were  used.  Assuming  that  the  third  version  of  the  coordination  function 
is  realistically  most  accurate,  there  is  a  tradeoff  between  the  amount  of  information  communicated  to  each 
ECAV  and  the  accuracy  of  the  team  cost  for  generating  a  phantom  track. 

Figure  10  below  illustrates  the  situation  once  the  variables  x$  and  a  as  well  as  f(x)  are  chosen  to  specify 
a  team-optimal  phantom  track  as  discussed  above.  The  bold  blue  line  represents  the  phantom  track  chosen 
optimally  via  the  coordination  functions  hi(f(x)),  and  the  yellow  lines  represent  the  resulting  range  of  valid 
initial  conditions  from  which  each  corresponding  ECAV  may  start  at  time  zero.  Time-dependent  flyable 
ranges  for  each  ECAV  are  not  shown  in  the  figure. 

Besides  using  coordination  functions,  several  other  methods  are  available  in  treating  the  decentralized 
cooperative  control  problem.  The  common  goal  is  to  use  the  minimum  amount  of  communication — or 
minimum  amount  of  the  team  state  information,  such  as  ECAV  positions — needed  to  achieve  the  desired 
level  of  performance  in  choosing  a  team-optimal  phantom  track. 
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Figure  10.  Decentralized  Cooperative  Control  Problem  for  3  radars  and  a  straight  phantom  track. 


X.  Conclusion 

This  paper  has  shown  that  for  mainlobe  deception,  assuming  the  ECAV  and  phantom  track  speeds  to  be 
constant  severely  limits  the  set  of  initial  conditions  yielding  complete  ECAV  trajectories  for  a  reasonable 
range  of  9  =  0  to  90  degrees.  Assuming  ranges  for  the  ECAV  and  phantom  track  speeds  and  including  a 
course  bound  to  account  for  the  ECAV  having  fixed  antennas  yielded  a  more  realistic  flyable  range  for  the 
ECAV.  Flyable  ranges  with  the  above  assumptions  were  thoroughly  investigated,  determined,  and  verified 
via  extensive  ODE  simulations  for  a  straight  phantom  track,  a  simple  circular  phantom  track  with  radar 
placed  at  its  center,  and  a  general  circular  phantom  track  with  radar  placed  arbitrarily.  The  approach 
developed  and  used  to  compute  ECAV  trajectories  and  bounds  given  a  straight  or  circular  phantom  track  is 
modular  for  n  radars  and  n  ECAVs.  The  ECAV  trajectory  solutions  require  specifying  one  ECAV  variable 
as  a  function  of  time — whether  it  be  speed,  course,  heading,  turn  rate,  speed  rate,  or  acceleration.  The 
general  bounds  developed  that  define  ECAV  valid  initial  conditions  and  flyable  ranges  were  shown  to  be 
independent  of  whether  the  phantom  track  is  straight  or  circular;  this  result  may  be  applicable  to  other 
types  of  phantom  tracks.  Based  on  these  general  bounds,  a  decentralized  cooperative  control  problem  was 
formulated  under  certain  assumptions  using  coordination  functions.  In  this  scheme,  communication  of  one 
coordination  function  per  ECAV — providing  its  individual  cost  for  generating  different  phantom  tracks — 
allowed  the  ECAV  team  to  determine  the  optimal  phantom  track  resulting  in  the  lowest  team  cost. 
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